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Background Symmetry restoration and configuration mixing in the spirit of the generator coordinate method based on energy 
density functionals have become widely used techniques in low-energy nuclear structure physics. Recently, it has been 
pointed out that these techniques are ill-defined for standard Skyrme functionals, and a regularization procedure has 
been proposed to remove the resulting spuriosities from such calculations. This procedure imposes an integer power of 
the density for the density dependent terms of the functional. At present, only dated parameterizations of the Skyrme 
interaction fulfill this condition. 

Purpose To construct a set of parameterizations of the Skyrme energy density functional for multi-reference energy density 
functional calculations with regularization using the state-of-the-art fitting protocols. 

Method The parameterizations were adjusted to reproduce ground state properties of a selected set of doubly magic nuclei 
and properties of nuclear matter. Subsequently, these parameter sets were validated against properties of spherical and 
deformed nuclei. 

Results Our parameter sets successfully reproduce the experimental binding energies and charge radii for a wide range of 
singly-magic nuclei. Compared to the widely used SLy5 and to the SIII parameterization that has integer powers of 
the density, a significant improvement of the reproduction of the data is observed. Similarly, a good description of the 
deformation properties at A ~ 80 was obtained. 

Conclusions We have constructed new Skyrme parameterizations with integer powers of the density and validated them 
against a broad set of experimental data for spherical and deformed nuclei. These parameterizations are tailor-made 
for regularized multi-reference energy density functional calculations and can be used to study correlations beyond the 
mean-field in atomic nuclei. 



PACS numbers: 21.60.Jz, 21.30.Fe, 21.10.Dr 

I. INTRODUCTION 

One of the most widely used designs of an effective 
nucleon-nucleon interaction for mean-field-based meth- 
ods [l| was introduced by Skyrme 0, [|[ as a combination 
of momentum-dependent two-body contact forces and a 
momentum-independent three-body contact force. Al- 
ready the first applications demonstrated the re- 
markable qualities of this interaction to describe many 
properties of nuclei throughout the chart of nuclei. How- 
ever, some drawbacks due to an insufficient flexibility of 
Skyrme's original ansatz became apparent. 

The early parameterizations of Skyrme's interaction 
led to two major problems. First, the simple contact 
three-body force does not allow for a realistic value of the 
incompressibility of symmetric infinite nuclear mat- 
ter. Typically, values of about 350 MeV were obtained, 
which is significantly larger than the empirical value of 
210 < K oo < 240 MeV 0-11 ■ Second, the same three- 
body force gives almost always rise to a spin-instability in 
infinite nuclear matter [10h12| | and finite nuclei [13j , ren- 
dering the calculation of excitations of unnatural parity 
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in RPA impossible jTI) . 

It turned out that both problems can be simultane- 
ously solved when replacing Skyrme's three-body force 
£3 8{r\ — r 2 )<5(ri — r 3 ) with a density-dependent two-body 
contact force 1 1 3 (1 + x 3 P a ) p%{(r 1 + r 2 )/2]S{r 1 - r 2 ), 

where P a and po( r ) are the spin-exchange operator and 
the isoscalar density, respectively. For a = 1 and X3 = 1, 
both are equivalent as long as time-reversal symmetry 
is conserved However, the so-called time-odd terms 
from the density-dependent two-body force have a differ- 
ent isospin structure than those of the three-body force, 
which removes the spin instability [10|. Therefore, all 
early parameterizations have since been used as two- 
body forces with a linear density dependence. In fact, 
the ambiguity around the three-body term was recog- 
nized from the beginning by the authors of the earliest 
fits of Skyrme's interaction 0, Hj], who pointed out that 
its three-body force "should not be considered as a real 
three-body force, but rather as convenient way of simu- 
lating the density dependence of an effective interaction" 

0- 

In a second step, reducing the exponent a to values 
between 1/6 and 1/3 allows also for a realistic compress- 
ibility 0, HH, [lH . The appearance of density dependen- 
cies of the form pfi (r) is also motivated through approx- 
imations to the G matrix of Brueckner-Goldstone the- 
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ory p7l - l20j . Up to now, all widely used parameteriza- 
tions of the Skyrme interaction have stuck to this simple 
form of density dependence, although several extensions 
were attempted over time [l| . The same form of density- 
dependent two-body contact force is also used to comple- 
ment the finite-range Gogny interaction [U [22j . 

Modifications of specific terms in the total energy have 
been made as well, hence abandoning the link to an un- 
derlying force Then, it is more appropriate to refer 
to a Skyrme energy density functional (EDF). 

There have been many adjustments of the parameters 
of Skyrme's interaction since the 1970s The range of 
data on which the parameters are fitted has been varied 
and extended, sometimes with choices dictated by spe- 
cific applications. Most fitting protocols, however, are 
designed to deliver multi-purpose parameterizations that 
can be used for applications as diverse as the descrip- 
tion of ground-state masses and density distributions, 
deformations, rotational bands, the response to external 
probes, fission, and reaction dynamics for nuclei all over 
the mass table and even the properties of neutron stars. 
One of these is the protocol by Chabanat et al. devel- 
oped in the 1990s [23|, [24|, which led to a significant im- 
provement of isospin properties by including pseudo-data 
for neutron matter. The resulting parameterizations, in 
particular SLy4, have been extensively tested and used 
for the description of many properties of atomic nuclei 
throughout the nuclear chart. 

The Skyrme interaction was designed for use in self- 
consistent methods, i.e. Hartree-Fock (HF), HF+BCS, 
Hartree-Fock-Bogoliubov (HFB), both in their static and 
time-dependent variants, and in RPA. From the 1990s, 
Skyrme's interaction was also frequently employed in ex- 
tensions of mean-field methods, such as the construction 
of a microscopic Bohr Hamiltonian [25H26j . exact projec- 
tion [271 128| and configuration mixi ng by the generator 
coordinate method (GCM) (2i| [3(], [32Tj35j . Two issues 
were apparent from the beginning, one related to the ad- 
justment of the parameters of the interaction and the 
other to its analytical form. 

Because the adjustment of the parameters is done 
within the mean-field approximation, the inclusion of 
beyond-mean-field correlations will often give rise to an 
overbinding of nuclei, in particular of those used during 
the fit. The extra binding, however, was always found 
to be within a few MeV and appears to saturate quickly 
when several collective modes and symmetry restorations 
are added consecutively [35[ . Here, we shall not consider 
this issue as overbinding constitutes a small smooth trend 
that at the present stage is smaller than other system- 
atic errors and/or uncertainties [3, H^, [36-38]. Instead, 
we will postpone the role of correlations on the outcome 
of a parameter fit to future work and concentrate on the 
setup of the functional itself 

As already mentioned, the Skyrme functional is 
adopted in extensions of the mean-field approach. Such a 
functional is a priori defined only for mean-field calcula- 
tions, i.e. for a single mean- field wave function, whereas 



beyond-mean-field calculations require to determine a 
matrix element between wave functions generated by two 
different mean fields. Unlike a formalism based on a 
Hamiltonian, the extension of a density functional from a 
single-reference (SR) definition to a multi-reference (MR) 
one is not an unambiguous procedure. In the early appli- 
cations of GCM using a Skyrme functional [3(|, the SIII 
and SIV parameterizations were indeed used as two- and 
three-body forces to calculate matrix elements between 
two different mean-field wave functions. Taking advan- 
tage of the generalized Wick theorem derived by Balian 
and Brezin [3l| , this amounts to replacing the mean-field 
densities that enter the energy density by so-called mixed 
densities. This scheme for the construction of the energy 
density was also followed in Ref. [32], where the param- 
eterization SIII was adopted as a density-dependent en- 
ergy functional to construct the energy kernel in a GCM 
calculation mixing mean-field states with different axial 
quadrupole deformation. This procedure was not altered 
until very recently in subsequent applications with more 
recent functionals of the Skyrme, Gogny and relativistic 
type [29|,|35|,[39|,[4(| that often include multiple symmetry 
restorations. 

However, such a generalization of the functional ig- 
nores several complications, in particular the fact that 
the mixed densities can become complex in MR calcu- 
lations and that, often, different functionals are chosen 
in the mean-field and pairing channels. Still, the sub- 
stitution of mean-field densities by mixed densities in 
the construction of the MR EDF was used with some 
success in many applications, despite its drawbacks that 
can in principle lead to unreliable results for an energy 
functional. In one way or the other, the problems of the 
functionals mentioned thus far are related to the breaking 
of the Pauli principle (4l| - |43l | . In the standard Skyrme 
EDF, this has many facets. First, the density depen- 
dence itself cannot be written in a completely antisym- 
metrized form. Second, it is customary to use different 
effective interactions for the particle-hole and particle- 
particle parts of the EDF. In addition, certain exchange 
terms of the Skyrme interaction are sometimes neglected 
or modified, and for the Coulomb exchange term ap- 
proximations are used. All of these are either motivated 
by phenomenology, or by computational reasons. For 
an overview, we refer to Ref. [l|. The standard den- 
sity dependence Pq(t) poses one additional problem. In 
all currently used prescriptions, the density entering the 
density-dependence might become complex. For non- 
integer values of a, the function Pq(t) then becomes a 
non-analytical function of pp that is multivalued and ex- 
hibits branch cuts j4l], SI, S3 • To resolve this particu- 
lar issue, some alternatives for the density dependence 
were formulated. Indeed, several studies have concen- 
trated on the most appropriate definition of the density 
dependence in MR calculations, primarily for symmetry 
restorations [H, SI, SB-SI]- The question, however, is 
not settled yet. 

The net result of these problems is that the off-diagonal 
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terms in the MR EDF can exhibit discontinuities or even 
divergences when v arying one of the collective coordi- 
nates. We refer to (ilUiiL Eij ] for an in-depth analysis 
of these issues but present the arguments for the lack 
of signs of their presence in the published GCM calcu- 
lations. First, the problems are especially critical for 
very light nuclei, but applications were often devoted to 
medium-mass and heavy ones. Second, the discretiza- 
tions commonly chosen for numerical reasons when set- 
ting up projection and GCM restrain the contamination 
of the energy with non-physical contributions to a very 
small scale. 

One possibility to avoid these problems altogether 
would be a return to a Skyrme-force-based Hamiltonian. 
This, however, will inevitably demand the systematic 
addition of higher-order terms in the Skyrme force, as 
within the standard form it is impossible to construct a 
parameterization that, at the same time, describes the 
empirical properties of nuclear matter, has no spin or 
other instabilities, and gives attractive pairing. By con- 
trast, within an energy functional framework a fair de- 
scription of nuclear matter and finite nuclei is achieved 
within the standard form. Thus, to keep the effective 
interaction simple, it appears to be preferable to work 
with a functional instead of a force. To enable their use 
in a MR framework, tools to by-pass the obstacles out- 
lined above by a regularization of the Skyrme functional 
have been designed recently [H, [4!|. They require, how- 
ever, that the functional dependence on the density has 
an integer power 44j . 

In this article, we construct Skyrme functionals that 
have the same density dependence as SIII and thereby 
are regularizable in the sense of Ref. [43]. The first pa- 
rameterizations of the Skyrme functional built about 40 
years ago [f| had all this property, but, since then, the fit- 
ting protocols have significantly evolved and these early 
parameterizations certainly have to be reconsidered. Our 
study is based on the protocol first used for the SLya; pa- 
rameterizations [23j, |24| that has proven to be efficient to 
construct functionals used successfully in a large number 
of applications. In this first study, we will restrict our- 
selves to the standard form of the Skyrme functional. 
The construction of a regularizable functional includ- 
ing higher-order density-dependent terms is underway 
[50l l5l| and will be reported elsewhere. However, we 
take the opportunity of the present study to include a 
new set of data in the fitting protocol, which are used to 
validate (or reject) the parameterizations. 

There is a major conceptual difference between the pa- 
rameterization of the Skyrme functional that we aim at 
and the ones by Kortelainen et al. [HI, [53|, who have 
recently adjusted new Skyrme parameterizations on a 
large set of data. The aim of Kortelainen et al. is to 
describe the nucleus in the spirit of the density func- 
tional theory [54| that is very successful in condensed 
matter physics. Staying on the computationally simple 
single-reference level, as much correlation energy as pos- 
sible is incorporated into the energy functional. Our aim 



is to construct a parameterization of the Skyrme EDF 
that will be used in beyond-mean- field calculations, i.e. 
where specific correlations are to be calculated explicitely 
in a multi-reference framework. Both views are com- 
plementary. The advantage of our approach is that it 
enables to calculate spectra and transition probabilities 
directly in the laboratory frame of reference and avoids 
the ambiguities related to approximate determinations of 
spectroscopic quantities, whereas its disadvantage is that 
already for standard observables high predictive power 
will require the time-consuming calculation of correla- 
tions beyond the mean field. In the following, we will 
call beyond mean-field method the method that we have 
already used in many applications and where mean-field 
wave functions generated by a constraint on a collective 
variable are projected on particle numbers and angular 
momentum and mixed by the GCM. 

The article is organized as follows. Section |TT] reviews 
the fitting protocol used here and its differences to the 
one used to construct the SLya; parameterizations in the 
past. In Sec. Mil we will test the parameterizations 
on a large set of typical observables for spherical and 
deformed nuclei, including masses, separation energies, 
charge radii, deformations, the fission barrier of 240 Pu, 
and the moment of inertia of a superdeformed rotational 
band in 1 Hg. Section HVl will summarize our findings. 



II. FITTING PROTOCOL 

A. The energy functional 

The standard density-dependent Skyrme interaction 
has the form [36] 

u(Rr) = t {l + x P a ) S(r) 

+±t 3 (l+Z3^)^(R)<5(r) 
+\t 1 {l+x 1 P a )[V 2 8{v)+8{v) k 2 ] 
+t 2 (l + x 2 P a )k' ■ 5{r) k 
+iW (<xi+<T 2 )-k'x£(r)k, (1) 

where we use the shorthand notation r = ri — r 2 and 
R = |(ri + r 2 ) for the relative distance and center- 

of-mass coordinates, respectively, where P a is the spin 
exchange operator, k = — |(Vi — V2) the relative mo- 
mentum operator acting to the right, and k' is the com- 
plex conjugate of k acting to the left, and po(R) is the 
isoscalar density. The Skyrme interaction (fTJ) contains in 
total 10 parameters io, tt, t 2 , t^, xq, Xi, x 2 , x$, Wo, and 
a to be adjusted to data. 

As it is customary, we only calculate the particle-hole 
part of the EDF from Eq. (p}. We keep, however, all 
terms in that channel, which is not always done For 
the special case of time-reversal invariance and spherical 
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symmetry this leads to 

^Skyrrnc = £ + + C^ P p t Ap t 

"* t=0,l 

+ Clptn + \c(Jl + c7 J Pt V ■ J t } , (2) 

where p, r, and J are the density, kinetic density, and 
spin-current vector density, respectively, and the index 
t labels isoscalar (t = 0) and isovector (t = 1) densi- 
ties. The definition of these densities and the relations 
between the coefficients Ct in Eq. ( El) a nd the parame- 
ters in Eq. ([I) can be found in Ref. [3g]. Note that the 
coefficients Cf[po] depend on the isoscalar density po{r), 
whereas all others are just numbers. In case of deformed 
nuclei and when breaking intrinsic time-reversal symme- 
try, there are additional terms in the Skyrme EDF for 
which we refer to Refs. [3g, |55| 

The total energy is given by the sum of the Skyrme 
EDF ([2]), the Coulomb energy, the kinetic energy, the 
center-of-mass correction and the pairing energy. As in 
our previous studies, we have chosen a density-dependent 
zero-range pairing interaction [l], [55|, [56[ , which leads to 
a functional of the form 



Vn 



£ fd*r 



q=p,n 



1 - 



Po(r) 



p q (r)p*Jr), (3) 



The switching density p c = 0.16 fin -3 is set to the em- 
pirical nuclear saturation density, such that the pairing 
interaction is most active on the surface of the nucleus. 
The pairing functional depends on the local pair density 
p g (r) 55] of protons and neutrons, labeled by q — p, n, 
and the isoscalar local density po(r). An energy cutoff 
of 5 MeV in the single-particle spectrum is taken above 
and below the Fermi energy [57|. The strength Vo will 
be adjusted separately for each parameterization of the 
Skyrme interaction. 

For most (if not all) Skyrme interactions constructed 
up to now, the Coulomb exchange energy has been re- 
placed by its Slater approximation [l| that amounts to a 
local energy density of the form ~ p^ 3 (r), i.e. a term 
depending on a non-integer power of the density. Like 
the standard density dependence in the Skyrme EDF 
~ Pq (r) with < a < 1, this term cannot be regularized 
with the currently available techniques [44j. For inter- 
actions that can be safely used in regularized MR EDF 
calculations, the Coulomb exchange energy has to be ei- 
ther treated exactly or to be omitted. For simplicity, we 
have chosen to neglect it in the mean-field channel in the 
present study since an exact treatment of the Coulomb 
exchange field makes all calculations much more time 
consuming. In addition, phenomenological arguments 
have also been brought forward that justify this course of 
action [58l - [60j . As usually done, the contribution of the 
Coulomb interaction to the pairing channel is neglected. 

For the center-of-mass correction, we employ the 
widely-used approximation where only the one-body 



term is considered [6l|] . However, the often neglected 
J 2 term in the Skyrme functional © is kept. The latter 
two choices correspond to the ones made for the param- 
eterization SLy5 of Chabanat et al. 24 1 . 



B. The protocol 

The first step of our fitting protocol is similar to the 
one used for the construction of the SLyx parameteri- 
zations (23l . |24| . During this step, we minimize a merit 
function which is a weighted sum of squared residuals: 



2 \ " 2 
X =2^,XA 



Xa 



1 



N A 



N A 

t—1 



AO, 



calc 



(4) 



The Oi are experimental data for finite nuclei and empir- 
ical values for nuclear matter and the AOi are tolerance 
parameters used to weight these data during the fit. Five 
categories of data are used: 

1. nuclear matter properties around the saturation 
point, 

2. neutron matter equation of state, 

3. binding energies of doubly-magic nuclei, 

4. charge radii, 

5. spin-orbit splittings of neutron and proton states. 
The nuclear matter properties that we have included are: 

• the saturation density p sa t = 0.16 fm" 3 with a tol- 
erance AOi = 0.003 fm -3 ; 

• the binding energy per nucleon E/A = —16 MeV 
with AO, = 0.3 MeV; 

• the symmetry energy a sym = 31 MeV with AOi = 
1 MeV; 

• the Thomas-Reiche-Kuhn sum rule enhancement 
factor k v = 0.25 with AO; = 0.15. 

Since the incompressibility of nuclear matter Koa cannot 
be adjusted to a realistic value with the restriction im- 
posed on a — 1 [23j ] . this quantity is not considered in 
our fitting protocol. 

The binding energies of six doubly-magic nuclei are 
included: 48 Ca, 132 Sn, and 208 Pb with tolerances of 
AO, = 0.2 MeV, 40 Ca and 100 Sn with AO, = 0.5 MeV, 
and 56 Ni with AO* = 0.75 MeV. We allow for larger 
AOi for N = Z nuclei as one always has difficulties to 
reproduce their binding energy at the mean-field level. 
However, the discrepancies cannot be simply related to 
the Wigner energy that cannot be described by mean- 
field calculations. Usually only 56 Ni turns out to be un- 
derbound, whereas 40 Ca and 100 Sn are overbound. The 
charge radii of 40 ' 48 Ca, 56 Ni, 132 Sn, 208 Pb have a toler- 
ance AOi = 0.02 fm, and the spin-orbit splittings of the 
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neutron 3p levels and the proton lh levels in 208 Pb have 
both a tolerance AOi = 0.2 MeV. 

For the neutron matter equation of state, Oi are the en- 
ergies per neutron for p 0.5 fm~ 3 predicted by Wiringa 
et al. [62] with the bare two-body UV14 potential and 
three-body UVII potential. The tolerance parameters 
are set to AOi — 0.2 x Oj. 

These data are used to determine a first set of val- 
ues of the Skyrme parameterization. The resulting EDF 
is then tested on several properties of finite nuclei that 
will be discussed in the following sections. Among these 
properties, the charge radii were strongly underestimated 
with the first set of weights that we have used. We have 
therefore chosen to relax the weights of nuclear matter 
properties, especially, the density at saturation and the 
constraints on neutron matter properties. After some 
attempts, this was sufficient to arrive to a satisfactory 
reproduction of charge radii. The weights that are given 
above are the final weights used in the fit. 

During the first attempts to fit our new parameter- 
izations, we encountered finite-size isospin instabilities 
that are characterized by a separation of protons and 
neutrons as examined in Ref. [631 ] . The instability ap- 
pears when the coupling constant C^ p — Jjti (| + X\) + 
^2 (5 + X2) in the Skyrme EDF ([2]) takes too large a 
value. To prevent such instabilities, we enforce a condi- 
tion on the coupling constant 



-c 



for C Ap ~> ^ Ap 



y l,calc 
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for C Ap alc < Ci 



Ap 
max' 



(5) 



where the empirical choice for the maximum value 
C A max = 25 MeV fm 5 has been found to lie safely within 
the stable zone. We have also checked that the param- 
eterizations do not lead to finite-size instabilities due to 
the Sj • Ast terms in the time-odd part of the Skyrme 
EDF [55j when setting the corresponding coupling con- 
stants to their Skyrme force value. 



III. RESULTS 



TABLE I: New parameter sets for the Skyrme energy func- 
tional with effective masses as indicated. 
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TABLE II: Saturation properties of nuclear matter as ob- 
tained with the new parameter sets. Values for SIII and SLy5 
are shown for comparison. 
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A. New parameter sets 

The fact that we do not constrain the compressibility of 
nuclear matter leaves some freedom in the choice of the 
effective mass, cf. the discussion in Ref. [23]. We have 
constructed four parameter sets corresponding to values 
of the isoscalar effective mass from 0.7 to 1.0 times 
the nucleon mass m. We will refer to these as SLylll.xa;, 
where xx is the value of m^/m. 

The coupling constants of these four parameterizations 
are listed in Table HI and the corresponding saturation 
properties of infinite homogeneous nuclear matter in Ta- 
ble [TTJ As expected, the value of is much too large. 
It increases with the effective mass [23| and there is no 
room to obtain a value close to the empirical value when 



FIG. 1: (Color online) Binding energy per nucleon E/A for 
symmetric nuclear matter (solid lines) and for pure neutron 
matter (dotted lines) as a function of nucleon density p. The 
filled squares indicate empirical values of the neutron matter 
equation of state by Wiringa et al. [63 |. 



imposing a = 1 without introducing additional terms in 
the Skyrme functional. To obtain a reasonable agreement 
between theory and experiment for charge radii has re- 
quired to relax the constraint on p sat leading to a value 
lower than the usual one of 0.160 fm -3 , but still larger 
than the one for SIII. 

The equation of state E/A of symmetric infinite mat- 
ter obtained with SLyIII.0.8 is compared with results for 
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TABLE IV: Pairing strength Vq for the parameterizations as 
indicated. The switching density is set to 0.16 fm -3 in all 
cases. 
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SLy5 and SIII in Fig. Q] As can be expected from the 
values for Kao, it is stiffer than the equation of state 
obtained with SLy5. 

In the same figure, we also compare the binding energy 
per neutron for pure neutron matter determined using 
SLy5, SIII and SLyIII.0.8 to ab-initio results obtained 
by Wiringa et al. [62]. On the scale of the plot, ob- 
vious differences between the parameterizations appear 
only at rather large densities p n > 0.12 fm -3 . At values 
below, the results obtained with the three parameter- 
izations cannot be easily distinguished, in spite of the 
fact that SIII was not fitted to this quantity. For larger 
densities, however, as expected, the inclusion of the neu- 
tron matter equation of state in the fitting protocol im- 
proves the results obtained with SLylll.xx with respect 
to those of SIII. For SLy5, the tolerance in the merit 
function, Eq. @ has been chosen much smaller than for 
the SLylll.xx, leading to a better reproduction of the 
equation of state. 

The residuals of binding energies and charge radii of 
doubly-magic nuclei are displayed in Fig. [2] and the corre- 
sponding values of % 2 are given in Tabic [TTT1 In both cases, 
the new parameterizations perform much better than SIII 
and SLy5, irrespective of the value of the effective mass. 
We have to recall, however, that SIII and SLy5 were fit- 
ted with different protocols, such that the comparison 
of the x 2 can on ly serve as a guideline for the relative 
performance of the parameterizations for these specific 
observables. It does not allow to judge their overall qual- 
ity. In particular, as discussed above, SLy5 gives a much 
better description of some key nuclear matter properties 
that cannot be adjusted with SLylll.xa:. 




4<V, 48,-, 56 xr 100 c 132 c 208™ 40„ 48^ 56 xr 132 c 208™ 
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FIG. 2: (Color online) Residuals between calculated and 
experimental values [M - [6r| . defined as 5 = (0 calc ' — 
O cxpt ')/|O cxpt |, for binding energies and charge radii of 
doubly-magic nuclei. 



TABLE V: Mean deviation with respect to the data for bind- 
ing energies and charge radii obtained for the parameteriza- 
tions as indicated. 





E dcv (MeV) 


E lel (%) 


^dcv (fm) 


r- r c e i (%) 


SLyIII.0.7 


1.97 


0.26 


0.018 


0.40 


SLyIII.0.8 


1.46 


0.21 


0.020 


0.46 


SLyIII.0.9 


1.09 


0.17 


0.023 


0.52 


SLyIII.1.0 


0.98 


0.15 


0.029 


0.65 


SLy5 


2.49 


0.31 


0.012 


0.29 


SIII 


1.88 


0.23 


0.051 


1.09 



erties of the interaction used in the particle-hole channel, 
we restrict ourselves to the surface pairing energy den- 
sity functional ^) that we have used in numerous past 
studies. 

The pairing strength Vq in Eq. ([2]) is fitted in 120 Sn on 
the neutron spectral pairing gap — £ pa iring,n/ Jd 3 r p n (r) 
[U [68|. In this expression, £ P airing,n is the pairing en- 
ergy of the neutrons and p n (r) the neutron pair density, 
respectively. The empirical value is determined by a five- 
point formula for the gap [13] and is equal to 1.393 MeV. 
The pairing strengths obtained for the four values of the 
effective mass are listed in Tab. IIVI They are very close 
to each other and do not scale significantly with the ef- 
fective mass. 



1. Adjustment of the pairing strength 

To compute spherical and deformed open-shell nuclei, 
pairing correlations need to be taken into account. The 
functional form and the adjustment of a pairing interac- 
tion is a problem that requires, in principle, a dedicated 
study of its own [67, 68] . Since our focus is on the prop- 



B. Spherical nuclei 

We start our validation of the SLylll.xx interactions by 
confronting their predictions with various experimental 
data for singly-magic nuclei. 
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FIG. 3: (Color online) Residuals of the binding energy E — £? e xpt. as a function of neutron number N for the Ca, Sn, Ni, and 
Pb isotopic chains obtained with the parameterizations as indicated. Experimental data are taken from Ref. [63 ] ■ 
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FIG. 4: (Color online) Same as in Fig. [3] but for N = 20, 50, 82, 126 isotonic chains as a function of proton number Z. 



1. Binding energies 

The differences between the calculated and the ex- 
perimental binding energies are shown in Figs. [3] and |4] 
for representative isotopic and isotonic chains of singly- 



magic nuclei. The agreement with the data is in general 
better for the SLyllLxx parameterizations than for SIII 
and SLy5. To quantify these energy differences, we have 



8 



defined two mean deviations: 



E, 



dcv 



1 N 



GXpt. I 



i=l 



E, 



rcl 



1 y> \E, - 



cxpt. I 



(6) 
(7) 



where N is the total number of singly-magic nuclei that 
have been calculated. Analogous quantities can be de- 
fined for charge radii. The values given in Table [V] con- 
firm that the agreement with data is improved by the 
SLylll. xa; parameterizations. Deviations for binding en- 
ergies -Edev decrease with increasing effective mass. 

Let us recall that our aim is to construct an interac- 
tion well suited for adding the correlations generated by 
symmetry restorations and configuration mixing calcula- 
tions. Therefore, the nuclei calculated at the mean- field 
level of approximation should be underbound, and that 
in such a manner that the difference between mean-field 
calculation and data is slightly larger for mid-shell nuclei 
than for doubly- magic ones [29j . 

It is clear that the SLylll. xx parameterizations with 
the largest values of m^/m leave nearly no room for the 
addition of correlation energies in the Sn and Pb chains. 
The increase of the effective mass washes out the shell 
effects in the mean-field results. At this point, SLylll 
TOp/m = 0.8 is the most promising parameterization, un- 
derbinding the energy of the Sn and Pb isotopes by what 
can be expected to be added from correlations. 



2. Charge radii 




112 116 120 124 128 132 
N 



80 82 84 86 
Z 



FIG. 5: (Color online) Charge radii for Ca, Sn, Pb isotopic 
chains and for TV = 50, 82, 126 isotonic chains obtained by 
SLylll, SLy5, and SIII. Experimental data shown by solid 
squares are taken from [6^ ]. 



The calculated and experimental charge radii are com- 
pared in Fig. [5] The charge radii are determined ac- 
cording to Ref. [HI, taking into account the internal 
charge distribution of both protons and neutrons and 
adding a correction for the electromagnetic spin-orbit ef- 
fect. The corresponding deviations, defined in Eqs. fl5J) 
and 0, are given in Table FVl The SLylll. xx parame- 
terizations clearly provide a better description of these 
data than SIII, which systematically underestimates the 
charge radii. SLy5, on the other hand, leads to even 
larger radii and therefore performs in general better than 
the SLylll. xx. As can be seen from Table fVl the devia- 
tions from the data i?dev and R re i decrease with decreas- 
ing effective mass. 

Again, we recall that correlations from fluctuations in 
the quadrupole degree of freedom consistently increase 
the charge radii of spherical nuclei [29| . Overall, the ob- 
served trend of the charge radii is well reproduced by the 
calculation. The deviations from the smooth trend ob- 
served in the data for the Pb and Ca isotopes and for 
the N = 82, and 126 isotones are not described by any 
of the parameterizations and seemingly require either the 
inclusion of explicit correlations, or higher-order terms in 
the EDF, cf. Ref. |l[ and references therein. 



3. Two-neutron separation energies 

The two-neutron separation energies are compared to 
the experimental data in Fig.|6]for the Ca, Ni, Sn, and Pb 
isotopic chains. All six parameterizations give similar re- 
sults for mid-shell nuclei. They tend to overestimate the 
characteristic jump at neutron magic numbers, which, 
however, would be reduced by dynamical quadrupole cor- 
relations [2{|. For Ca and Ni isotopes, our values do not 
reproduce the slope of the experimental data for mid- 
shell nuclei. For the Sn and Pb isotopes, the agreement 
with the data is improved by SLylll. xx with respect to 
SIII and SLy5. 



4. Single-particle energies 

Up to now, our analysis of the Skyrme parameteriza- 
tions has been limited to data for which the comparison 
between theory and experiment is model independent. 
This is no longer the case for single-particle energies, for 
which there exist several conflicting definitions that often 
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FIG. 6: (Color online) Two-neutron separation energies for 
Ca, Ni, Sn, and Pb isotopes calculated with the parameteri- 
zations as indicated. 



mass. The SLy5 results are sometimes very different. In 
all cases, the reproduction of the experimental gaps is 
rather poor. A more detailed analysis would require to 
compute directly one-nucleon separation energies, includ- 
ing correlations beyond the mean-field which are known 
to give a sizable contribution to the two-nucleon sepa- 
ration energies to and from doubly- magic nuclei p9l [73[ . 
We present below in Sec. IIII CMl results for self-consistent 
calculation of binding energies of a few very heavy odd-^4 
well-deformed nuclei, for which correlations beyond the 
mean field can be expected to play a lesser role. 



C. Deformed nuclei 

In addition to the properties of singly-magic nuclei, we 
also validate the performance of the new parameteriza- 
tions for deformation and rotational properties of selected 
key nuclei. 



1. Deformation energy curves 

We start by studying nuclei that have been experimen- 
tally identified to be either deformed or have states of 
different deformation coexisting at low energy. In Fig. [31 
the deformation energy curves of 24 Mg, 74 Kr, 80 Zr, 100 Zr, 
and 186 Pb are plotted as a function of the dimensionless 
axial quadrupole deformation 



even do not correspond to observables [70j • Here, we use 
here the eigenvalues of the single-particle Hamiltonian. 
They provide a lowest-order approximation to separation 
energies, which should be corrected for polarization ef- 
fects [H and the coupling to vibrations, cf. Refs. [7ll. l72l| . 
In Fig. single-particle energies are compared to one- 
nucleon separation energies to or from doubly-magic nu- 
clei. 

For 40 Ca and 132 Sn, the spectra obtained with the 
SLylll.ra parameterizations are very similar to those of 
SIII. For 208 Pb, there are several differences, in particular 
concerning the position of high-j levels. 

A rule of thumb predicts that a higher effective mass 
gives a more compressed spectrum. We are in a good po- 
sition to check this rule. The SLyllLzx have been con- 
structed using exactly the same protocol but correspond 
to four values of the effective mass. They are a mod- 
ern version of SIII but share with it many similarities. 
Looking to Fig. a higher effective mass corresponds in- 
deed to a more compressed spectrum. However, a change 
in the effective mass does not correspond to a simple 
rescaling of the single-particle spectra. For neutron holes 
in 208 Pb or neutron particles in 132 Sn, the relative dis- 
tances between levels hardly vary at all. Also, this rule 
of thumb is already not valid anymore for a change in 
the fitting protocol, as exemplified by SIII. The differ- 
ences between the single-particle spectra obtained with 
SIII and the SLyllLccx cannot be due to the effective 
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16tt 3R 2 A 
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(8) 



where R — 1.2 A 1 ' 3 fm. Experimentally, a prolate defor- 
mation for their ground state is well established for 24 Mg 
and 100 Zr. For 80 Zr, spectroscopic data suggest that the 
excited states of the ground state rotational band have 
a large quadrupole deformation with a 02 value around 
0.4. The sparse available data, however, do not rule out 
that the ground state of 80 Zr has a complicated structure 
that involves a large mixing of different deformations. 

The only parameterization that gives rise to a pro- 
nounced prolate minimum for 24 Mg, 80 Zr and 100 Zr is 
SIII. For both SLylll.xa; shown, the ground state of 80 Zr 
is spherical with a prolate minimum at a slightly higher 
energy, nearly degenerate with a very shallow oblate min- 
imum, whereas the ground state of 100 Zr has a large pro- 
late deformation, with an oblate minimum at smaller 
|/32 1 excited by around 1 MeV. For SLy5, the ground 
state of 80 Zr is spherical with a prolate minimum ex- 
cited by about 4 MeV. For 100 Zr, this parameterization 
gives nearly degenerate prolate and oblate minima. Be- 
fore drawing conclusions on how well these topographies 
are compatible with experimental data, one has to es- 
timate how the correlations that we plan to introduce 
explicitly in future applications might change the simple 
picture of energy curves. Rodriguez and Egido [74j have 
calculated the energy surface of 80 Zr including triaxial 
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FIG. 7: Single-particle energies obtained as eigenvalues of the mean-field Hamiltonian in spherical calculations for 40 Ca (left), 
132 Sn (center), and 208 Pb (right) for the parameterizations as indicated and compared with experimental data determined as 
one-nucleon separation energies to or from the doubly-magic nucleus. 




FIG. 8: (Color online) Potential energy curves along quadrupole deformation fa for 24 Mg, 74 Kr, 80 Zr, 100 Zr, and 186 Pb obtained 
by constraint HFB calculations. 



quadrupole deformations using the Gogny force. They 
have found several spherical, axial and triaxial minima. 
Before projection, the axial part of their energy surface is 
similar to the one obtained here with SLy5. For this nu- 
cleus, however, projection on angular momentum alters 
the topography of the energy surface, leading after con- 
figuration mixing to a ground state with a predominant 
component at a large quadrupole axial deformation. 

A beyond-mean-field study of the neutron-deficient Kr 
isotopes using the Skyrme parameterization SLy6 has 
been published in Ref. |34j. After projection and mixing, 
the relative energy of prolate and oblate states leads to 
excitation spectra in disagreement with the experimen- 
tal data. The energy curve obtained for 74 Kr with SLy5 
resembles the one of SLy6 presented, such that it can be 



expected that SLy5 would also give similar results after 
configuration mixing. By contrast, the prolate minimum 
obtained with SIII and SLylll.xr seems more realistic in 
view of the experimental data. Finally, the deformation 
energy curve of 186 Pb is alike for all parameterizations, 
displaying a spherical ground state and a prolate and 
oblate minimum within less than 1 MeV excitation en- 
ergy each. 

It is remarkable that SIII and the SLylll.xa; parame- 
terizations give a much more realistic description of the 
energy curves in the A w 80 region than the SLycc param- 
eterizations. This difference, however, cannot be traced 
back directly to the linear density dependence, as some 
other Skyrme parameterizations with non-integer expo- 
nents a of the density dependence give an energy curve 
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FIG. 9: (Color online) Fission barrier along quadrupole defor- 
mation /?2 for 240 Pu obtained by different Skyrme parameter 
sets. 



for 80 Zr that is much closer to the one SIII than that of 
SLy5 (z|. 

The value of the effective mass has a clear effect on 
the variation of energy with deformation. Comparing 
the curves obtained with SLyIII.0.8 and SLyIII.1.0, one 
can see that a higher effective mass results in a flatter 
behavior of the deformation energy curves. 

Overall, the SLylll.xa; parameterizations provide en- 
couraging results for the deformation properties at the 
mean-field level. The following examples, however, will 
illustrate some limitations of these parameterizations. 



2. Fission barrier of 240 Pu 

In Fig. HI the fission barrier of 240 Pu is presented as a 
function of quadrupole deformation. For all parameter- 
izations, triaxiality was taken into account in the calcu- 
lation of the first barrier and octupole deformations for 
the second barrier. 

In all cases shown in the figure, the excitation energy of 
the fission isomer overestimates the experimental value, 
for which two conflicting values of ~ 2.8 MeV [76| and 



2.25±0.20 MeV 77] can be found in the literature. In the 



same way, the energies of the inner and outer fission bar- 
riers overestimate the experimental values of 6.05 MeV 
and 5.15 MeV respectivel y [781] . For SIII this deficiency 
has been known for long [79|]. Also, the results obtained 
with SLy5 are less realistic than those obtained with the 
SLy4 and SLy6 parameterizations discussed in Ref. [80l |. 
However, one must take into account that the calcula- 
tions performed in Ref. [8(| and here are not fully equiv- 
alent: the pairing strength is not the same and particle 
number projection was performed in [80| and is not here. 
In that paper, it was shown that at the mean-field level, 
the energy of the fission isomer is close to the experimen- 
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FIG. 10: (Color online) Dynamical moment of inertia as a 
function of angular frequency hui for the superdeformed rota- 
tional band in 194 Hg obtained with the parameterizations as 
indicated. 



tal value with SLy6, whereas SLy4 gives better agree- 
ment when beyond-mcan-ficld correlations are taken into 
account. All parameterizations shown in Fig. [S] give en- 
ergies for the fission isomer much larger than SLy4 when 
used with standard pairing, and beyond-mean-field cor- 
relations cannot be expected to be large enough to obtain 
agreement with the data for any of them. 

The differences in barrier height for SIII and SLylll.ia; 
seen in Fig. [9] cannot be correlated with the value of the 
surface energy coefficient of these parameterizations. The 
values for SLy5 (a sym = 18.5 MeV) and SIII (18.6 MeV) 
are very similar, whereas those for SLyIII.0.8 (19.5 MeV) 
and SLyIII.1.0 (19.4 MeV) are significantly larger. The 
value of the isoscalar effective mass, and thereby the av- 
erage level density at the Fermi energy, does not play 
a crucial role either. However, the similarity of the en- 
ergy curves obtained with SIII and all SLylll.ia; hints at 
an insufficiency of a simple linear density dependence to 
describe large deformation. 



3. Superdeformed rotational band in 194 Hg 

The next test of the new parameterizations concerns 
the description of superdeformed rotational bands (SD). 
These bands are well described all over the nuclear chart 
by self-consistent mean-field calculations and represent 
one of the most impressive successes of these approaches 
in the 1990's. The SD bands in the Hg region are of 
specific interest as the gradual increase of the dynamical 
moments of inertia J^ 1 



J 



(2) _ 



d(jj 



1 d£_ 



(9) 



as a function of rotational frequency fiw results from 
the gradual disappearance of pairing correlations and the 
alignment of the intruder orbitals. For further details we 
refer to our recent detailed analysis of the various contri- 
butions of the EDF to in Ref. [Uj]. The dynamical 
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FIG. 11: (Color online) Low-lying one-quasiparticle states in 
251 Cf and 249 Bk. Experimental data taken from [H[I1- 



moments of inertia for the ground SD band of 194 Hg are 
presented in Fig. [TQl For SIII and SLyllLzx, the peak in 
the J^> appears at too low an hw and, overall, the de- 
scription of the experimental data is less satisfactory than 
that of other Skyrme parameterizations such as SLy4 or 
SkM*. It should be noted that the currently used pairing 
strength Vo is rather low in comparison with the typical 
values of Vo = 1250 MeV fm 3 that was determined in SD 
bands. An increase of the pairing strength, however, will 
have very little influence on the location of the peak in 
3^. 



4- Single-particle levels in deformed transactinide nuclei 

Previous studies of odd-mass transactinides [82| have 
put into evidence some major drawbacks in the spec- 
tra obtained with the current Skyrme parameterization 
s. We have tested the parameterizations that we have 
constructed in this work on two nuclei 251 Cf and 249 Bk, 
for which very detailed data are available and which have 
been studied in Ref . 82] . The same method as in Ref. [82[ 
has been used. Each state results from a self-consistent 
calculation of a one-quasiparticle excitation on an even- 
even HFB vacuum. In this way, the polarization effect 
due to the quasiparticle excitation and the terms in the 
Skyrme EDF depending on time-odd densities are taken 
into account self-consistently. The results are shown in 



Fig.QTJ The 251 Cf spectra exhibit the expected effect of 
the effective mass: the spectrum is becoming more dense 
when the effective mass is increased. Note, however, that 
the compression of the spectrum is not uniform and that 
the changes do not correspond to a simple scaling propor- 
tional to the ratios of effective masses as it is sometimes 
assumed [85j]. Moreover, the order of the levels can be dif- 
ferent when comparing the parameterizations. The non- 
trivial effective mass dependence is still more apparent 
for the spectrum of 249 Bk, where the first excited state 
is lower in energy for the lowest values of the effective 
mass and does not have the same quantum numbers for 
all the parameterizations. Although the obtained spec- 
tra depend on the parameterizations, none of the SLy.xx 
corrects the main drawbacks of previous EDF parameter- 
izations, i.e. the misplacement of some levels that may be 
connected with specific spherical single-particle orbitals. 



5. Particle number symmetry restored deformation energy 
surfaces 

The main motivation of the present study was to con- 
struct a Skyrme functional that can be used in regular- 
ized MR EDF calculations. In Fig. [T2J we show as an ex- 
ample of such a calculation the particle-number restored 
deformation energy surfaces of 24 Mg without and with 
regularization in one sextant of the ,8-7 plane. The cal- 
culations were performed as described in Ref. [4!| with 
two differences. The first one is the use of SLyIII.0.8, 
and the second one is the use of the extension of the reg- 
ularization scheme to trilinear terms in the same particle 
species as is required by this parameterization. We use 
Fomenko's prescription [49] with 19 discretization points 
for the gauge-space integrals. At small deformation, the 
difference between the regularized and non-regularized 
energy surfaces is quite dramatic. Without regulariza- 
tion, the absolute minimum is located in a region where 
the spurious contribution to the EDF is particular large. 
It is only with the regularization that one finds the usual 
topography of the energy surface with a prolate axial 
minimum. 

The nature and size of problems with spurious contri- 
butions to the MR EDF depend strongly on the param- 
eterization of the functional. The presence of terms that 
are trilinear in the same particle species in SLyIII.0.8 
makes the deformation and discretization dependence of 
the spurious energies much more violent than what is 
found for the SIII parameterization used in the regular- 
ized calculations Ref. [49j]. Also, there are no evident 
problems in the (non-regularized) particle-number pro- 
jected energy surfaces of 24 Mg obtained with SLy4 and 
presented in Ref. [35[ • There, we encountered obvious ir- 
regularities only when projecting simultaneously on par- 
ticle number and angular momenta J > 0. 

A detailed discussion of the regularization that will 
also address its application to angular-momentum pro- 
jection and general configuration mixing will be given 
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FIG. 12: (Color online) Particle-number restored deformation 
energy surface of 24 Mg in the /3-7 plane constructed with the 
SLyIII.08 MR EDF without (panel a) and with regulariza- 
tion (panel b) of the MR EDF. Black filled circles indicate 
the location of the absolute minimum of each energy surface, 
and grey circles the location of secondary minima. The dots 
indicate the calculated points. 

elsewhere [86j |. 

IV. SUMMARY 

The present study is a part of our program to con- 
struct an effective interaction of high spectroscopic qual- 
ity for mean-field and beyond-mean-field calculations. In 
this first step, we have constructed a regularizable (in 



the sense of (43[) EDF within the standard form of the 
Skyrme EDF. This requires that the power a of the den- 
sity dependence takes an integer value. The simple form 
of the non-momentum dependent trilinear terms used 
here has known deficiencies. It forbids to obtain a value 
for the incompressibility compatible with the empirical 
value. We have shown also the problems encountered in 
the description of charge radii, fission barriers heights and 
moments of inertia in SD bands in the A » 190 region. 
However, the protocol that we have developed leads to 
a significantly improved description of shape coexistence 
in the A w 80 region. 

The four variants with different isoscalar effective 
mass will enable studies on how the correlation energy 
in beyond-mean-field methods depends on the effective 
mass. However, the mean-field results presented here 
show a clear preference for mj/ra = 0.8. 

Even with their deficiencies, the present parameteriza- 
tions will allow us to analyze and benchmark the perfor- 
mance of the regularization. Work in that direction is 
underway fS6j | 

Higher-order terms {i.e. at least trilinear terms with 
derivatives) are clearly necessary to remove the defi- 
ciencies of the SLyllLix parameterizations pinpointed 
here and to improve the predictive power of regulariz- 
able Skyrme- ty pe f unctionals. Work in that direction is 
also underway [50l l5l| . Alternative (non-Skyrme-type) 
regularizable forms of the density dependence might be 
considered as well, cf. for example the form proposed in 
Ref. [I?]]- The moment the form of a sufficiently flexible 
functional that is safely usable in MR EDF calculations 
has been established, fits should be performed on the 
level of MR EDF. 
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